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Abstract. In this paper we consider a system of equations that describes a class 
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stable mechanisms of singularity formation (see e.g. M. P. Brenner et al. 1999, 
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prevails over the diffusion and a non-self-similar shock wave results. Our main 
result identifies a class of initial data for which the blow-up behaviour is of the 
former, self-similar type. The blow-up profile is characterized as belonging to a 
subset of stationary solutions of the associated ordinary differential equation. 
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1. Introduction 



We consider the parabolic-elliptic system 

nt = div{eVn + n'Vcj}} 
A(j) = n 



m 



m 



n X M+, 
n X M+, 

dn X IR+, 
dn X IR+, 

n, 
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n{x, 0) = nQ{x) 
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where = -Bi(O) = {a; € M'*: < 1}, c? > 2, and i7 is the outer normal vector from 
the boundary dQ. Here > is a constant parameter. The initial condition uq is 
chosen in L^(f2), radially symmetric, and such that 




and 




(6) 
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Equations C3l~© define a problem for the unknown mass density n and potential 4>. 
Mass is conserved by the no-flux condition and therefore © implies 

n{x,t)dx — J no{x)dx^l. (7) 
a SI 

Problem CQl-© is a model for the evolution of a cluster of particles under 
gravitational interaction and Brownian motion (see [S] and the references therein). 
Here n represents the mass density, (jj the gravitational potential, and O a rescaled 
temperature characterizing the Brownian motion. This model also appears in the 
study of evolution of polytropic stars, by considering the evolution of self-interacting 
clusters of particles under frictional and fluctuating forces (22j. Finally, problem 1^- 
(O also arises in the study of the motion of bacteria by chemotaxis as a simplification 
(see of the Keller-Segel model 122 EH El IH- Here the variables n and cf) represent 
the density of bacteria and the concentration of the chemo-attractant. 

We view the problem Q)-® as an evolution equation in n, since by equations 
OJ the function (j> is readily recovered from the solution n. It is known |H] that problem 
(C3)-||njl has a unique local solution if Uq S L^(17), which satisfies n £ L°°(f2 x (e,T)) 
for some T > and for every e > 0. We restrict ourselves to the analysis of radially 
symmetric solutions and write n{r,t) := n{x,t) with r = |a;| € [0, 1]. 

Since we are interested in the question when and how 1^-® generates 
singularities, we define: 

T* = sup{ r > I Problem iP)-® has a solution n e L°°(fi x (e, t]) }. 
If T* < oo, then we say that blow-up occurs for Ql-®, in which case 

lim sup n{r,t) = oo. (8) 

*^^* [0,1] 

Various sufhcient conditions for blow-up are known O^KTHS]. 

For d ~ 3, Herrero et al. |19[ OU] were the first to study the behaviour of the 
solution close to blow-up, using matched asymptotic expansions. Later Brenner et 
al. studied the problem for 2 < d < 10. They used a numerical approach to 
describe solutions and proved existence and linear stability of similarity profiles. Note 
however that no proof of convergence or characterization of blow-up in terms of initial 
data were given in these references. The principal types of blow-up described in 
[iniEnilini are: 

(a) A solution n{r,t) consists of an imploding smoothed shock wave which moves 
towards the origin. As t T*, the bulk of such a wave is concentrated at 
distances 0((T* - t)'^/'^) from the origin, has a width 0((T* ~ and at 

its peak it reaches a height of order 0((T* — t)"2(d-i)/d^^ rpj^jg ^^^^ blow-up 
has the property of concentration of mass at the origin at the blow-up time, i.e. 



lim 

r— >o 



^lim, y n{y,t)y'^ ^ dy 





= C > 0. (9) 



This situation is depicted in Figure ^ (left), 
(b) A solution n(r, t) has a self-similar blow-up of the form 

(T* - t)n (jjy/iT* - t)e, tj - ^{tj) as t ~* T* . (10) 

Note that this implies that n satisfies © with C = 0. Therefore no concentration 
of mass at the origin occurs at the blow-up time. This blow-up behaviour is 
depicted in Figure Q (right). 




The results of this paper are two-fold. First, we demonstrate rigorously that the 
self-similar blow-up structure (|10|l is an attractor for the system secondly, we 

identify an explicit class of initial data that converges to a self-similar solution of this 
type. Let us elaborate on this. 

Let no = no{r) be such that 

Xdr'^mir) < ||no||Li(B,,(o)) for r £ (0, 1), (11) 

Q{no)r + noiMr > 0, (r'^(0o)r)r = r'^no in (0, 1), and M'^) = 0, (12) 

where Xd is the measure of the unit ball in R'^. Suppose also that 9 < l/(4dxd), 
implying that the solution n = n{r,t) of Q-® blows up at finite time T* > and 
at the point r = Finally, assume that the two functions 

||?^o||li(b,-(o)) and 2(d - 2)9r* + ^'^^'^'^^'^'^^ exactly once in [0, 1]. (13) 
Our main result f Theorem 12. If) shows that if (|1HI . (|12|l . and l|13|l hold, then n satisfies 

and moreover has a structure near blow up given by 
n,(r,i) = (T* -t)-i* 



Ve(r*-i); 

where the function 5* is one of a class of solutions of a steady-state problem; a class 
that includes the functions 

^^J,{r^):=id-2)-^^±4\^^ and ^*{r^):^l for ,7 > 0. 
(rf - 2 + ^r]^y 

In particular the initial state ng = l/xd and O < l/(4c?Xd) satisfies the conditions 
above (Corollarv l2.2|l . If we relax assumption (|13|) but assume instead that n satisfies 
the growth condition 

n(0,t) < M(r* -t)-i for te(o,r*), 

for some constant M > 0, then n has the same structure of blow-up given above 
(Theorem I2.3|) . The hypotheses on the initial data (fTT|l . ({T^ . and ifT^ are more 
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natural in the context of a transformed problem we introduce in the next section. 
Note however that (no)^ < in [0, 1] implies assumption 

This paper is organized as follows. In section[21 we put the problem in terms of a 
new variable, thus transforming the system into a single PDE, and then state 

our results in terms of this new formulation. In section |31 we discuss some non-self- 
similar blow-up patterns related to case (a). Sections 01 El and El provide the tools for 
the proofs of Theorems 12.11 and 12 . 31 and the arguments are wrapped up in Section d 
A rather technical dervation of a Lyapunov function is placed in Appendix and in 
Appendix 151 we derive some linear stability results. 



2. Precise statements of main results 



For radial solutions, the average density function b{r,t) ^Hj is defined by 

r 

b{r,t) :^ f niy,t)y'-' dy, (14) 



This variable turns out to be convenient in the analysis of this system. Note that it has 
the same scale invariance as n(r, t), but that solutions are smoother when expressed in 
terms of b. For example, if for some fixed t > the density n(r, t) is a delta function at 
the origin with unit mass, then b{r,t) — r^"^. Let D = (0, 1) and set Dt = D x (0,T) 
for some T > 0. Transformation (|14|l puts system (Q-© in the form 

bt = Xd& (^brr + -7-^'-^ + ^rbbr + b^ in Dt (15) 

br{0,t) = 0, 6(l,i) = l, for te[0,r), (16) 

6(0,0 = 60(0 for reD. (17) 

Here we have redefined t := —t. Regarding the initial condition, wc assume 

60 e C^(D), and -(bo)r + feo > for r e D, (18) 
d 

where the second condition is equivalent to no > in _D. Note that the conservation 
of the mass O is represented by b{l,t) = 1 for i G [0,T). As was done for problem 
(C3)-||ni) we define T > to be the maximal time of existence for the average density 
b{r,t). If T* < 00 in ©, then 

lim sup b{r, t) — 00, 

where T = T* /xd- Using (fnji . we deduce b(r,t) < for r e f, t > 0; this implies 
single point blow-up for 6(r, t) at the point r ~ {). To characterize the asymptotic 
behaviour near blow-up of the solution b{r,t) of problem p5(l - (|18|l . we study the 
solutions of the associated boundary-value problem 

Vrm + + - -jJl'Pv + <^ - V = 0, for 77 > 0, ^^^^ 

^(0) > 1 </7„(0) = 0. 

If 6 is a solution of p5|) - (|18|l which blows up at time T > and at the point r = 0, 
then we will show that it has the asymptotic form given by 

-1, 



5,(r,i) = (r-t)-V 
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Equation ifT^ has multiple solutions for 2 < d < 10 HOI El- We classify them by 
counting the number of times they cross the singular solution (psiv) ■— 'id/rf'. For 
that purpose, we introduce the set 

Sk = {</3: is a solution of (|19|l that has k intersections with tps}- 

We shall see that iSi is the relevant subset of solutions of p9|l for the characterization 
of the type of blow-up considered in this paper. Numerical evidence QUI suggests that 
iSi contains only two elements: 

^*{r,)^l and ^1(77) := ^ for > 0. (20) 

(d 2 + 2 J 

For the initial condition, we assume 

(6o)r < for reD, (21) 



and 



Xde[{ba)rr + ^^ibo)A +^rbo{bo)r + bl>0 for r e D. (22) 



We will show that this implies br < in Dt and &t > in Dt- In terms of no 
assumption (|21() becomes (lll|l and assumption (|22|l becomes (|12|) . 

Theorem 2.1 Let d > 2 and bo satisfy l^21\) and Igjj)) . Let b{r,t) be the corresponding 
solution of problem i)_?5|) - i)-?<j^) that blows up at r — and att — T. If 



< Oi := l/(4(ixd) o.nd ba{r) intersects T ^Lpi{r / yxd^T) once (23) 
then 

6(0,t) < Mi(r-t)-i for te{Q,T) (24) 
with Ml :— 2d/{d— 2). Moreover, T < Mi/6o(0), and there exists f £ Si such that 

^lim (T - t)b {v^XdQ{T-t)) = ^(77) (25) 

uniformly on compact sets jryj < C for every C > 0. 

We remark that there exists a family of 60 satisfying the conditions H18|l . H21(l . and 
H22|l . given by bQ{r) = Ki + K^jir'^ -f -^^3) with positive constants Ki that satisfy 
Ki + K^lil + i^g) = 1 and 9 < fCz/Srf^Xd- Conditions 10, HH), and igll) are also 
satisfied for &o = 1- Note that condition (|23() of Theorem 12.11 can be generalized 
by changing (pi for other solution (p of p9(l . Since these solutions are only known 
numerically, the counterpart of Mi and Oi cannot be given explicitly. The next 
corollary applies this result to 60 = 1- 

Corollary 2.2 Let d > 2, bo = 1, and Q < Qi. Then b{r,t), the corresponding 
solution of problem il_?.5|) - |77[ ), blows up at r = and at some time t ~ T < Mi; 
m,oreover {24^ holds and there exists ip (£ Si satisfying l^25\) . 

Numerical simulations 10 suggest that for an open set of initial data the convergence 
in (|25|l holds for ip — ipi. This self-similar behaviour may be seen roughly in Figure ^ 
(right), by imagining n{r, t) replaced by b{r, t) (since n and b scale similarly). InlBlwe 
show that ipi is linearly stable (using the result in ^Hl) and also that ip* is linearly 
unstable. 

For more general initial data we have the following result. 
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Theorem 2.3 Let d > 2 and let fog satisfy l^21\) and {S^. Assume that b{r,t), the 
corresponding solution of problem I^15p - fl8p . blows up at r = and at t — T. If b 
satisfies the growth condition 

b{0,t) < M{T -t)-^ for <e(0,T) (26) 

with M > 0, then there exists € 5*1 such that the convergence 125\) holds. 

We now briefly discuss the structure of the proofs of these theorems. Following 
the scale invariance, we set 

r = log(^), ^- (^^e(T-t))V2 ' and Bi.^r) ^ (T - t)bir,t). 
The rectangle Dt transforms into 

n = t) \ T > 0, < ri < e{T) } where ^(r) := {xdQT)-'^/^e''^^ . 
The initial-boundary problem p5l) - (|18|l now becomes 

Br + B+ ^rjB,j = Br,r, + '^-^Br, + \vBB,-, + £2 in n, (27) 

B^(0,r) = 0, B{1{t),t) = e-^T forTe]R+, (28) 

B(77, 0) = Bo{fi) T&o (r^ixd^TY'^) for r; G n(0), (29) 

where 11(0) = (0,£(0)). Note that a solution of l|19|) is a time-independent solution of 
H27|l -p9 |) . Therefore the study of the blow-up behaviour of b{r,t) is reduced to the 
analysis of the large time behaviour of solutions B{ri, r) of (|27|I - H29|) . and in particular 
stabilization towards solutions (p of (|19|l . The proof of Theorem 12.31 consists of two 
parts. In Sectional we first prove that lo <Z Si, where 

uj = {(j)£ L°°(R+) : 3Tj oo such that 

B{-,Tj) ^ (p{-) as Tj ^ oo uniformly on compact subsets of R^} (30) 

is the w-limit set we introduce for (|27|I - H29|I . The proof uses the observation that 
equation (|27|l . without the convection term i 77 _Bi?,j, is the backward self-similar 
equation for the parabolic semilinear equation 

bt = ANb + b^, (31) 

where A jv denotes the Laplacian in and N = d-\- 2 ^1 ^1 . We use the methods 
for the analysis of this self-similar equation to prove Theorem 12.31 However, due to 
the presence of the convection term, a different Lyapunov functional is necessary. This 
functional is constructed using the method of Zelenyak PU], which yields a Lyapunov 
functional in implicit form. In section (BJ we use intersection comparison arguments 
based on the ideas of Matano to prove that the w-limit set l|5n|) is a singleton. 
With a result on intersection with (ps this completes the proof of Theorem 12. 31 

Note that Theorem l2.3l is similar to a result for the supercritical case {N > 6) for 
equation H31|) , where two different kinds of self-similar blow-up behaviour may coexist 

m 

Finally to obtain Theorem 12.11 and Corollary 12.21 we use Theorem 12.31 and 
comparison ideas from Samarskii et al. 26 , Chapter IV] . 
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3. Discussion on non self-similar blow-up patterns 

In this section we discuss a family of blow-up patterns which appears when we refine 
the asymptotic expansion for the profile (p — (p* = 1. This situation is closely related 
to the blow-up behaviour of H31|) with < 6. If a solution b of (|31|l with < 6 blows 
up at a: = and t — T, then 



hin (T - t)b{r]y/T ~ t, t) = 1 

uniformly on compact sets |?7| < C for arbitrary C > Moreover it have 

been shown (see for instance [23123) that a refined description of blow-up gives the 
existence of two possible types of behaviour: either 

^lim (T - t)b (vViT-t)\\og{T^t)\, t) = ^1 (77) (32) 

uniformly on compact sets \ri\ < C, with C > arbitrary; or 

hrn (T - t)b (tj(T - ty^^"',tj = (^^(ry) for some m > 2, (33) 

uniformly on compact sets \ri\ < C, with C > arbitrary. Here the family {ifi}i>i 
is known explicitly. For problem H15|l - (|18|l . it was shown [20] for d ~ 3 that there 
exists a refined asymptotics for (p* = 1. Extending the argument to all d > 2, these 
asymptotics suggest a convergence given by either 



,y^(r-t)|iog(r-t)|('^-2)/^ ?j=p 



^lim(r-t)6(77V(r-t)|log(r-t)|('^-2)/'i,i) =p,{r^) (34) 
or 

^lim (T - t)b (ri{T - i)^+2,„^^i, ^ = ^„(ry) (35) 

for some m > 2. An implicit formula for the family {(^m}m>i is given in |1(JI 
equation (43)]. The type of convergence in 77 towards these profiles is an open problem. 
In H35|) , we can take formally the limit m — > cx) and find a non-trivial scaling, 

^lim (T - t)&(?7(r - ty/'', t) = poo{v). (36) 

Note that this limit cannot be taken for the semilinear equation where H33|l holds. The 
convergence H36f) represents the convection-dominant behaviour of (|15|) - (|18|l . which in 
terms of the density n = n(r, t) describes an imploding wave moving towards the origin, 
as shown in Figure ^ (left). The function 1^00 is discontinuous (cf. (3.16)]), 

— — for 7] > C 



for < C, 



where 2C"^ is the mass accumulated in the origin, which can be chosen arbitrarily. 
In this type of blow-up was studied using matched asymptotic expansions. There 
it was suggested that this behaviour is stable and moreover it was expected that there 
exist initial data such that (|36|l holds uniformly in rj on compact subsets away from 
the shock. A result of this type was proved in |12[ Theorem 3] for a related equation. 
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4-1. Estimates 

In this section we develop some estimates for problem (|15|l - (|17f) . which in turn will 
imply bounds for the self-similar problem H27f) - (|29|l . 

Lemma 4.1 If bo satisfies ^S\) then 

-hr + b>0 m Dt- (37) 
a 

Proof. The solution n of problem satisfies the relation 

1 r 

n^—\-X + b] in Dt-- (38) 
Xd a 

Since no > in Z?, an application of the maximum principle to problem JQl-lini) shows 
that n > in Dt- - Using this and H38() the result follows. ■ 

To prove the following results, we proceed as in 113! where similar estimates were 
found for the semilinear parabolic equation (j31|l . 

Lemma 4.2 //&o satisfies \21\) then 

br.{r,t)<0 m Dt- (39) 

Proof. Set w{r,t) :— r'^^^br{r,t). Differentiating (|15|) . we find 

Wt ~ XdQ (wrr ~~''^r] - '^^^'^r ^ ( ^ + d'"^'' ' ^' ^"^^^ 



r 

Assume for the moment a stronger assumption on the initial data 

(6o)r(r)<0 for rG(0, 1) and (6o)rr(0) < 0. (41) 

This gives ■w{Q,r) = r'^^^br{0,r) < 0. Under H41() the function 6 = 1 is a sub-solution 
for (|15|I - H18() . but not a solution; by Hopf's Lemma, w{l, t) — br{t, 1) < for all t > 0, 
so that w < on Dt, hence 6^ < on Dt- To finish the proof, we note that by the 
strong maximum principle, if bf) satisfies (|21ll . then for each ti G (0,T) condition ()41|l 
holds for the function b{r,ti). This proves the result. ■ 

Lemma 4.3 // bo satisfies 1^21]) and assuming that blow up occurs at time T > , 
then 

b{0,t)>{T-ty' for te[0,T), (42) 

Proof. Since the maximum of b in D is attained at r = (by Lemma |4.2I) . we have 
brr{0, t) < 0. It follows from (|15() that 6f (0, t) < 6^(0, t). Integrating this inequality on 
(0, T) gives the result. ■ 

Lemma 4.4 If bo satisfies H^) then bt > for all t e (0,r). 

Proof. Condition (|22|l implies that bo is a subsolution for I|15I17II : therefore 6(r, e) > 
b{r,0) for small e > 0. By the comparison principle we find b{r,t -I- e) > b{r,t) for 
t e (0, T - e). It follows that 5* > on Dt- ■ 
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The next lemma gives a bound on \bj.\ in D^- 
Lemma 4.5 Let bo satisfy if^) and if^j . T/ien 

Xde52(r,i) < ^6(0,i)3 /or (r,t)ei^T. (43) 
Proof. Since bt > and 6^ < in Dt, we multiply equation (|15fl by and obtain 

r 



= IxMblir, t) - bl{Q, t)] + ^b^r, t) - ha\o, t). 
Since &^(0,t) = we obtain the desired inequahty. ■ 

To conclude this section we translate the properties of solutions derived above 
into estimates for problem (|27|I " H29|I . From hypothesis 1)26(1 and noting that 6 > 1 and 
6r < in DT^ we have the a priori bound 

< 5(7/, r) < M for (r;, t) G H. (44) 
Combining this with (|43|) and 139|) . we obtain 

< -5^(77, t) < M for (?7,r)en, (45) 
where M depends on M . Finally from 142|l . we get 

l<B(0,r) for Te(0,^(T)). (46) 

T/ie steady state equation ^19p 
We begin by recalling problem H19|l : 

V»?r; + ~~~'-Pn + - ^Wri + (fi^ - If = for > 0, (47) 

^(0) > 1, ^,,(0) = 0. (48) 

Condition (|48|l is required, since 5(0, t) > 1 for all r > 0. Equation 147|l has three 
special solutions: 

fsiv) = -^, <i5*(?7) = l, and (p^{r]) = for 77 > 0. 

Note that (^55 satisfies 

+ ^v{'^s),i = and = ((^35)^^ + ^^^^(^5),, + ^77(^5(^35)^ + (y's)^- (49) 

For bounded non-constant solutions we have the following theorem ^1 1201 ■ 

Theorem 4.6 Let 2 < d < 10. There exists a countable set of solutions {ipk}keN of 
j-^7| )- |7^ such that ipk{0) > 1 and ipk{0) — > 00 as k ^ 00, Moreover ipk intersects the 
singular solution ips k times and has the asymptotic behaviour (pk{rj)rj'^ = Const{k) > 
0. 
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The proof is based on the equation for G{ri) :— r]'^ip{r]), 

hm < oo, Urn -nG^ir]) = 0. (51) 

Note that ips corresponds to G{'i]) = 2d. 

It was formally argued in 10 that for each integer k > 2 and 2 < d < 10 the set 

Sk = f solution of H47|l - (|48|l with k intersections with ^ps} 

is a singleton and that for d > 2 the set Si contains only two elements. More precisely, 
iSi consists of the functions (p* and ^pi given in (|2()|l . If we relax condition (|48|l to 
</5(0) > 0, we conjecture that there is at least one other solution in Si. For d = 3 this 
was shown numerically by Brenner et al., who found a solution Lp\ of l|47|l such that 
(/3*(0) < 1 and {^*i)^{0) = 0, which intersects (ps once ^ Figure 14]. 

5. Convergence 

In this section we prove the following convergence theorem. 

Theorem 5.1 Let conditions and hold. Let B{rj,T) he a uniformly hounded 
global solution of \21^ - \29\) . Then for every sequence t„ — > cx) there exists a 
subsequence such that B{r],T^) converges to a solution ip of j^'^ ^ f?^ - The 
convergence is uniform on every compact subset o/[0, oo). 

Proof. Define B'^{ri,T) := B{ri,(T + r). We will first show that for any unbounded 
sequence {uj} there exists a subsequence (renamed {%}) such that B"^ converges to 
a solution ip of (|4 7(1 -148 (1 uniformly in compact subsets of M+ x M. Without loss of 
generality we assume that the sequence {uj} is increasing. 

Let N G N. We take i large enough such that the rectangle Q2N — {{v^''') G 
M.'^-.O < ri < 2N, \t\ < 2N} lies in the domain of B"' . The function B(^,t) = 
B"' (I'f I, t) is a solution of 

Br = Ad+2B - • VB + ^(e • VB)B + B^-B 
on the cylinder given by 

TsAT = {(e,T):R'^+2 X M: |e| < 2N, |r| < 27V}, 

and \B{^,t)\ is uniformly bounded in r2Ar by (|^ . 

By Schauder's interior estimates all partial derivatives of B can be uniformly 
bounded on the subcylinder C T2n- Consequently B"' , Br' , B^' , and S^'^ are 
uniformly Lipschitz on Qn C Q2N- By Arzela-Ascoli, there is a subsequence {nj}f^ 
and a function B such that B"', i?"', B^S and _B^j^ converge to B, Br, B,j, and i?,,,,, 
uniformly on Qjv- 

Repeating the construction for all N and taking a diagonal subsequence, we can 
conclude that 

S". ^ B, B^^' ^ Br, ^ B^, and B"^ ^ B^^, (52) 



Self-similar blow-up 11 

uniformly in every compact subset in M+ x M. Clearly B satisfies (|27|) and estimates 
H44|l and (|45|l . Finally, it remains to prove that B is independent of r. This implies 
that i? is a solution of H19|l . since B{0, r) > 1 for all r > 0, and the result follows. 
Claim. The function B is independent of r. 

To prove this we construct a non- explicit Lyapunov functional in the spirit of 
Galaktionov and Zelenyak 1301 . 

1. Non-explicit Lyapunov functional. We seek a Lyapunov function of the 

form 

£(r) 

E{t)^ J <^{'q,B{'q,T),B^{n,T))dj^, 



where €(t) = {xdQT)-^/'^e'' /"^ and $ = $(r/, v,w) is a function to be determined. In 
Appendix A we show that such a Lyapunov function exists; more precisely, we show 
that a function p = p{r],v,w) exists such that 

e(r) 

d f 

— E{t) = - j p{7j, B{r,, T),B.,{r,, r)) {Br)^V, r) drj 



+ <i>»S.ir' + lKr)^{i{T),B{l{r),r),Br,{l{r),r)). (53) 

To identify the relevant domain of the functions $ and p, we note that by estimates 
(gl)) and (gni the solution B satisfies {r/, B{r/,T), B,j{r/,T)) e TZ, with 

'R.^nn{0<v<M,0<-w<M}, (54) 

where 7^ {?] > 0, w > 0, w < 0} U {?? = 0, w > 0, w = 0}. 

The functions p and $ are continuous in TZ \ {rj — fj , v > 1} with fj > defined 
later and they satisfy 

l^d+lg-Cor,^ < p{rj, V, W) < r^<l+l^-{d-2)rf/4d f^^. ^ -j^^ (55) 

Co 

with Co = Co (A/) > (Lemma ESJ, and 

|$(?7, w)| < c^^d+ie-(^-2)')V4rf i;, w) e (56) 

for some positive constants Ci{M) > fLemma IA.6|I . 

2. Proof of the claim. An integration over the interval {a,b) of H53|l gives 

b 1{t) 

p(ry, B(ry, r), i3^(r;, T))i?2(^, r) drydr i?(a) - i?(6) + V(a, 6) (57) 



a 

where 



6 

1 



^{a,b):= / -£(T)$(f(r),B(£(T),r),B,(€(T),r))dr + 



6 

+ / i?.(£(T),T) 



a 



B,(£(r),r) 

p(£(r),B(^(T),r),s)ds 







dr. (58) 



Since 5,(£(r),r) = -S(£(T),r) - i^(T)B^(^(T), r), 

i?.(£(T),T) = -Te-^ ~ \br{l,T{l - e^)). 
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Applying (|37|l at r = 1 gives 16^(11^(1 — e"^))! < d and consequently is uniformly 
bounded as r — s- cx). Employing this bound on B^- and the estimates (|55|l and (|56|l we 
find 

lim {supV'(a,6)} = 0. (59) 

<^'^°° b>a 

By (|52|) . we have that there exists a sequence nj oo such that B"^ [t], t) converges 
to B uniformly in compact subsets of (K+)^. For any fixed N we will prove for a 
subsequence satisfying lim (n^+i — rij) ~ oo that 

hm / p{rj,B"^{rj,T),B:;^{T^,T))iB!^^)'{7j,r)drjdT^0, (60) 

rtj — >oo J ' 
Qn 

where we recall that Qat = {(»y, t): R^: < r/ < iV, |t| < N}. Since p is bounded from 
below on bounded subsets of TZ, it then follows that 

[ Bldr]dT= lim [ iB^^f{r],T) drjdr ^ 0, 

J nj^oo J 

Qn Qn 

proving the claim. For all j sufficiently large, 

N < (xdeT)"i/2e^("^-^) and n^+i - Uj > 2N. 
Consequently using H57|l . we find 



N N 



-N 



p{7i B"^ (77, r), (77, r))(B;'^ )2(r;, r) d,7dr 

< / / p{r^,B^^^{^,T),B;^{rj,T)){B^^)'{f^,T)dvdi 



-N 



< 

ni~N 



p{Tl, B{rj, r), 5,(77, r))(S0'(?7, r) ^77^7 



< E{nj - N)- E{nj+i - N) + 7/^(71^ - N, n^+i - N). 
Hence applying l(5^ . we discover 

p(77, 5"^ (77, r), 5;'^ (77, t)){B"^^ t) dvdT < limsup[i?(77, - N) - E{n,+i - N)]. 



Qn 

Next we divide the expression £'(77,^ — N) — E{nj+i — N) into three integrals, choosing 
K arbitrarily large: 

E{nj ~N)~ E{nj+i - iV) = 

K 

- I [$(77, B"^ (77, -N),B;^ {in, -N)) - $(77, S"^ (77, -iV), B^^^ (77, ~N)] drj (61) 


+ y $(77,S"^+i(77,-7V),5;^+n^,-^))'^'7 (62) 

K 
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T-l/2g 2" 



+ / <i>(7?,i3"^(r/,-7V),i3;^(?7,-7V))(i77. (63) 



K 



Integral H61|l tends to zero as j ^ oo. In fact by the continuity of $ in the second and 
third argument we obtain pointwise convergence and by the bounds H56|l on $, we 
apply the Dominated Convergence Theorem to conclude. Expressions and 
can be made arbitrarily small since they can be bounded by 

oo 



K 



where C is a positive constant, and K can be chosen arbitrary large. Thus we have 
proved Hfi()|l . concluding the proof of the Theorem. ■ 



6. Comparison results 

6.1. Comparison with the singular solution Lps 

This section closely follows 1 . From section l?!^ we recall that solutions </? of H47l) - (|48|l 
are classified by their intersections with ips- In this section we study the intersections 
of solutions B of (|77|l - lf^ with ips- Our results are closely related to the ones found 
in PP, where equation (|31|l was studied. 

We first see that for 8 < l/{2dxd) a solution B of H27|l -p9 |l intersects the singular 
solution ips at least once in 11(0) since 

(^s(O) = oo > B(0, 0), and ps (iXdQT)-^") < B (^(xdOT)-'^^ , o) = T. 

On the other hand, for 8 > l/{2dxd) it can also be shown that B intersects (ps at least 
once in 11(0). Assuming the contrary, suppose that B{-,0) < ps{-) in 11(0). By the 
maximum principle, we obtain B < ps in 11. Therefore in the limit r — > oo, thanks to 
the Theorem 15.11 and since B{0, r) > 1 for all r > 0, we find a solution ip of (|19|l such 
that (p < ips- However we can show that every bounded non zero solution p of (|19|l 
has to cross (ps- This is equivalent to proving that there exists no solution G of (jSOII - 
H51|l such that G{ri) < 2d for 77 > 0. To check this, we assume that such a solution 
exists; we examine two cases. Suppose that for some rj*, we have Gri{r]*) — and 
G{r]*) < 2d. By (I50f) . G has a strict minimum at 77*, which contradicts the boundary 
condition (|51l) . On the other hand if, G{ri) is increasing for all 77 > 0, then for large 
77, equation H50|l implies that Grjrj > 0, which also contradicts (|5HI . 

We conclude that there exists rji S 11(0) such that -6(771,0) = (ps{T]i) and 
B{r],0) < ips{ri) for rj < 771. 

Lemma 6.1 Under the assumptions \21]) and f^jj)) . there exists a continuously dijfe- 
rentiable function 7yi(T) with domain [0,oo) such that ?7i(0) = r]i and B{7]i{t),t) = 
fsivii'^)) for all T > 0. 

Proof. Define H{r],T) := B{ri,T) — ipsiv)- We first claim that H,Hjj, and Hr do not 
vanish simultaneously. Using Lemma |4 . 41 and the strong maximum principle we find 

bt^iT-t)-^ fBr + B+^7^BA >0 in Dt. (64) 
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Suppose there exists a point in 11 where H,^ = = H = 0. Then iJ^ = imphes 
Bt — 0, and condition H,j — combined with H = gives 

B + ^jjBr, = in n, 

using H49|l . This imphes that 6t = at some point of Dt, a contradiction with Ht)4|) . 
Secondly, we claim that H^i 7^ at any point (fi,f) G 11 where H(fi,f) = and 
moreover H{ri, f) < in a left neighborhood of fj. A proof of this can be done as in 
Moreover, from the proof, we find -ff^(?7,f) > 0. 

Now we prove that 7J^(?7i,0) > 0. This follows from the equation satisfied by 
H{tj, 0). To the left of we find 

^^^^(^7,0) + ^H„(77,0) + ^T^H^iTj,0){B{7^,0)+^s) (65) 

+ j^vH{v,0){B{^,0)+^s)v>0. (66) 

Since (-8(77,0) + ifs)rj < and H{rji,Q) — 0, we can apply Hopf's Lemma to obtain 
that Hr,{r]i,0) > 0. Finally, to conclude the proof of the lemma, we use the implicit 
function theorem as in ■ 

Define the set Hi = {(77, t) | < 77 < 771 (r) } and the function 

mir) = . gyp|^ g (^^^ {xdOT)-'^^]:H{s, 0) > for s e [771, 77]}. 

Since H{rii,0) = and iJ^(?7i,0) > 0, the above supremum is finite. Define the set 

n2 = {(77, t) I 7/1 (r) < 77 < 772(t)}. 
Let F{t) = H{f]2{T),T). By definition of 772, F(0) > 0. Also, 

^F{t) = HMt),t) + ^772 (r)i/, (7,2 (t),t). 

Using lfi4|l . we have dle"^ F (t)] / dr > 0. An integration yields F{t) > for t > 0. 

As was done in pp, applying the maximum principle, using Lemma l6.ll and noting 
that H{r]2{T),T) > for r > 0, we can prove the following lemma and its corollary. 

Lemma 6.2 The function H{ri,T) ~ B(r],T) — ips{ri) satisfies H < m Hi and H > 
in 112- 

Corollary 6.3 Assume the conditions in Lemma \6. 1[ For each N > there is tn > 
such that for r > rjv, ^(77, r) intersect (psiv) most once m 77 G (0, TV). 

6.2. Intersection comparison 

In this section we derive comparison results, which will be used to prove that uj, the 
limit set (jSOJ, is a singleton. 

We start by considering the following linear equation with inhomogeneous 
boundary conditions: 

Vt — Vrr H — — Vr + for 0<r<l, ri<i< T2, 

r 

Vr{Q, t) = for Ti < i < T2, 

v{l,t) = h{t) forTi<i<T2; 
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where Ti , T2 are positive constants and 

aei°°([0,l] X (Ti,T2)), h e C\{Ti,T2)), (68) 
are given functions. Moreover we assume 

h{t) > for Ti < t < T2. (69) 
The zero number functional of (|67|l is defined by 

z[v{;t)] = #{re[0,l]:v{r,t)^0}, (70) 

and the foUowing lemma provides some properties of this zero number functional. 

Lemma 6.4 ([24 ) Let v = v{r,t) be a nontrivial classical solution of \61\j and 
assume that \6b]) and hold. Then the following properties hold true: 

(i) z[v{-,t)] < 00 for any Ti < t < T2; 
(a) z[v(-,t)] is nonincreasing in time; 

(Hi) ifv{ro,tQ) — Vr{ro,tQ) = for some tq G [0, 1] and to > Ti, then z[v{-,t)] drops 
strictly at t = to, that is, z[v(-,ti)] > z[w(-,<2)] for any Ti < <i < tg < t2 < T2. 

From this lemma we deduce a property of intersection between a solution ip of (|19|l 
and a solution B of |[77|l -ll^. 

Lemma 6.5 Let B be a bounded solution of | |i^7| )~ il£9)) and let ip be a solution of 
Denote Z{t) = fj={r G [0, ^(r)]: _B(?7, t) = (^(77)}. Then the following properties hold 
true: 

(i) Z{t) < 00 for any t > t*; 
(a) Z{t) is nonincreasing in time; 

(Hi) if B{r]o,To) = ipivo) and B,^{rio,To) = <y£'r,(??o) for tq > n, and 770 < ^(t) then 
Z{ti) > Z{t2) for any n < tq < T2. 

Proof. Writing V = U - b, where U(r, t) ^ {T - t)-^ip{r/{xdQ{T - t))^^'^), we 

have 

Vt = Vrr + ( + -u] Vr + (-X + b + u) V for < r < 1 , < t < T, 



(71) 

K(0,0 = 0, V{l,t)^U{l,t)~b{l,t) forO<i<T. 
Let Ti < T2 < T. For the variable F(r, t) = exp ^ yU{y, t) dyj y(r, t), we find 

Vt = Vrr + ^^^Vr + A(r, t)V for < r < 1, Ti < t < Tj, 

r 

Vr{0,t) = 0, for Ti<t<T2, 

V{l,t)^{U{l,t)-l)exp^^lyU{y,t)dy'^ for Ti < i < Ta, 
where 

Air, t)='fr + b + U+l[ yUtiv, t) dy - ±,r^U^ - (t/ + rUr) - '-±^U. 
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Note that A e i°°([0, 1] x (Ti, Ta)) since b, 6^, U, Ut, Ur £ i°°([0, 1] x (Ti, Tz)). If we 
show that V(l,t) does not change sign for t > to, then setting 7\ = and using 
Lemma 16.41 we have proved the lemma. 

We claim that there exists to such that Ut{l,t) does not change sign for t > to. 
By definition of V, this implies that there exists to > such that V{l,t) does not 
change sign for t > to- 

Since Ut{r,t) = (T - t)-'^{ff(p),J{2ri), if r = 1 and t > t\ then 

Ut{l,t)^{T-t)-^^(Tf^)^ for t>t\ and ii>rj*{t*), (72) 

where r?*(t*) {xd^iT — t*))^^/"^ . From Lemma A.l], we know that for a given 
a G (0,4(i), any solution ip of H47|) satisfying 

ri^ip{ri) a as r/ cx), (73) 

is such that there exists fjo — fjo{o-) so that the sign of {ri'^ip)^^ does not change on 
[770,00). Using H72|) . this implies that there exists to = ^0(^0 ) such that the claim 
holds. ■ 



7. Proofs of main results 

We start by proving that the a;-limit set of problem H27|I ~ H29() is a singleton. 

Theorem 7.1 Assume the hypotheses of Theorem \2.Sl Then the setuj defined in ]30(l 
is a singleton. 

Proof. For this proof we extend a solution B of (|77|l - lf^ to aU (M+)^ by setting 
B{r],T) = e^'^T for (77, r) G (M+)^ \ H. We also define the weight function p*{r]) ~ 
e-')V4 foj. ^ > 0. 

The hypothesis (|26|l implies that B is uniformly bounded; Theorem l5. II therefore 
states that uj is non-empty, and that each ip £ uj is a. solution of (|47|l ~148 () . 

We claim that for each ip £ oj there exists r* > such that B{0,t) — ip{0) never 
changes sign in [t*,oo). By contradiction, we assume that there exists a sequence rj,, 
such that Tfe 00, and B{0,Tk) — (p{0). Since i?^(0,Ti;) — fij{0) — 0, by Lemma IH^ 
the function Z{t) has to decrease at least by one. However this cannot happen an 
infinite number of times. This proves the claim. 

Suppose now that lu is not a singleton. Since the w-limit set is connected, 
closed, and non empty, it contains an infinite number of elements. We select three 
different elements tpi,ip2, pa in the cj-limit set. Since these functions are different and 
each solves ()19|l . we may assume that </3i(0) < ^'2(0) < PsiO)- By the claim above, 
B{0,t) — ip2{0) never changes sign in [t*,oo). This contradicts the fact that ipi and 
(P3 are elements of w; it follows that u; is a singleton. ■ 

We now conclude the proof of Theorems 12 . 31 and 12 . II and Corollarv l2.2l 

Proof of Theorem 12.31 By the previous theorem w is a singleton, say {B}. 
From CoroUarv 16.31 we find that for every > there exists a r^f > such that 
the solution 5(77, t) intersects fsiv) most once in r] £ [0,N] for each r > tn- 
This implies that in the limit t 00, B intersects ips at most once, concluding the 
proof. ■ 
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Proof of Theorem [Zn Since b and Ui{r,t) = (T - t)" Vi(r/(x<i6(r - t))^/^) 
are solutions of p5|l with the same blow up time, V — b — Ui satisfies equation (|71|) . 
Using that Ui{r,t) = 2d/{{d - 2)iT - t) + / {2xdQ)), we find 

V{l,t) = {l-Ui{l,t))>0 if e < l/(4dxd) for any t<T. 

The functions Ui with 6 necessarily intersect exactly once for all t, since non- 
intersection implies that the solutions must have different times of blow-up |26l p. 
271]. It follows that b{0,0) < Ui{0), and one finds {T - t)b(0,t) < 2d/{d- 2). An 
application of Theorem 12.31 proves the theorem. ■ 

Proof of Corollary If 6o = 1 and 6 < l/(2(d -I- 2)xd), we know 

from j?! Theorem 2] that the corresponding solution b blows up. Now assuming 
6 < l/(4dxd) < l/(2(d + 2)xd), we can apply Theorem lO to conclude. ■ 



A. Appendix: The Lyapunov functional 

In this appendix we construct the Lyapunov functional E satisfying l|53(l with the 
suitable properties of p and $ to prove Theorem 15.11 We start with a formal 
construction of the functional. This requires solving a first-order equation for p after 
which $ can be expressed in terms of p. Finally, we explain how to use smooth 
approximations of $ to obtain a rigorous derivation of l|5r 



A.l. Formal derivation of a Lyapunov functional 

Assume that $ and p are regular. To find such functions satisfying H53(l . we compute 

£(r) i{T) 

■^EiT)= J $,B,rfr/+ J 'P^Br^dr)+^-^^i{T),B{£{T),T),B^i£iT),T)). (A.l) 



Wherever possible we omit the arguments of <& and p for clarity. Integrating by parts 
the second integral in ljA.l|l becomes 

e(r) 1(t) 

r) 



j <^>niBrr, drj = - J [$^^ + <^vwB,, + ^niwB,„]Br dr] + ^^S^jg^' 



Defining 



ff \ V , ^ ,2 

/ ( 7y, w, w) = IV w H — r]viv -\- v — v, 

?/ 2d 



equation H27|) takes the form = B^^ + f{rj, B, B,^), by which equation ljA.l|) becomes 

i{T) 

d f 

-^E{t) ^ / {[<^v-^,jro-<^vwB,,+^njwf]Br-<^ww{Br)^}dr] 


+ <f>.B.|f ) + ^M^e^r),Biiir),r),B,ie{r),T)). 

Now if functions p = p{ri,v,w) > and $ = exist that satisfy the system 

of equations 

-^v + ^-r/u, + w<^^^ = pf and ^^ow = P, (A.2) 
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then E has the form of a Lyapunov functional with a contribution on the boundary, 
i.e. 

eir) 

Mr) 



d f 

■^E{t)=- J p{rj,B,B„){Brfdii-^<^>^BrC- 





+ -\imr),B{£{r),T),B,{e{r),r)). (A.3) 



2 

Therefore we may obtain this formula by solving system ljA.2|l . which we do by 
transforming it to a first-order equation for p, 

wpv + p,, - fpw = fwP- (A. 4) 

If we supplement a given solution p of this equation with the function $ given by 



$(77, v,w) ^ J {w - s)p{t], V, s)ds- J p{r], p, 0)f{r], p, 0) dp, (A.5) 



then the pair (p, $) solves (|A.2|I . In order to find the pair (p, $) we therefore only 
need to solve equation ljA.4|) . 



A. 2. The first- order equation for p 



We solve equation ljA.4p by the method of characteristics. Characteristic curves of 
equation (|A.4p are curves x = {ri,v,w) in M.^, which we consider parametrised by 77, 
along which 

d d / A ^\ 

—v — w and —w = —f. (A. 6) 

d-q dr] ^ ^ 

If a curve x(77) = {r],v^{ri),w^{ri)) satisfies these equations, then equation (|A.4|I 
reduces to 

^p(x(r;)) = /^(x(r7))p(x(ry)). (A.7) 

In order to solve the system of ODE's (|A.6p and (jA.7|l . we select a vector (770, fo, wq) G 
E+ xM^ and define (/)(^) = 0(f ; 779, vq,wo) to be the solution of the initial value problem 

0" + /(^, 0') = 0, with (l)\^=no=vo and (f>'\^=^g = wq, {A.8) 
where ' = If the curve x passes through {r]Q,vo,wo), i.e. if x(7;o) = {tlOiVojWo), 
then this curve can be identified with 770, wq, wq), since x(77) = {r] , [r]) , (rj)) 
where 

v^iv) = Hv;Vo,vo,wo) and ^^(77) = ^'(77; 770, wq, ^o)- (A. 9) 

Since = — ^ + ^777;, we may integrate (|A.7|I to find 

n 

p{r],v,w) ^ p{r]o,vo,wo)exp\ [ ^ 2 ^ ^^^'^ 

p{Vo,vo,wo)'^e~^'^'+^o/iexph J ^v^Odd- (A.IO) 



Vo 
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To prove Theorem 15 .11 we need to define p in the sellZ ClTl given by H54|l . 
7^ = {?? > 0, u > 0, w < 0} U {77 = 0, w > 0, w = 0} 

n = nf^ {0 < V < M, < -w < m}. 

We do so in the following way: for each (ry, v, w) € 7?., we define p(?7, f , ly) by following 
the characteristic curve through (ri,v,w) to a reference point (770, uq, wq) for which 
p(?yo, I'D, Wo) is fixed by choice; the value of p(r], v, w) is then given by iTOO)! . To select 
an appropriate set of reference points, we study some of the properties of solutions (j) 
of iTOt . since they define the characteristic curves. 

It follows from standard ODE theory that solutions of (|A.8|I are locally smooth 
and continuous under changes of {r/o,VQ,wo). In general, however, we cannot extend 
these solutions to the whole of M"*"; in fact, for each {ri,v,w) € TZ, there may exist 
< ^1 < and/or ^2 > V such that 

<p{^i; r/,v,w) — 00 and/or ^(^2; ?7, w) = — 00. 

Partly because of this difficulty, we choose to only use forward solutions of l|A.8(l to 
define the characteristic curves. The next result details the behaviour of a forward 
solution (j) of 1A.8|I . 

Lemma A.l Let {ri,v,w) G TZ, and let 0(^) = (j){^;ri,v,w) be the solution of fA.^) . 
For ^ > rj, exactly one of the following three alternatives holds: 

(i) 0=1 or (f) = 0; 

(ii) there exists rj* > rj such that 4>{rf) — and 4>{Ct < for ^ > r]*', 
(Hi) (piO ~* as ^ ~> 00 and there exists a constant C > such that (f>{£,)(,'^ —> C as 
^ 00. 

Proof. See ^| p. 95]. The proof is based on results from [20]. ■ 

Since we need to define p with the appropriate estimates, we introduce a 
parameter fj in the following lemma. 

Lemma A. 2 There exists fj > such that for every rji > fj any solution (j) of ifvl.iS)) 
with 4>{rji) — 1 and 0'(?7i) < satisfies 

(p'im) < -1 for all T]2 > ?7i with (/-(yya) E [0, 1/2]. (A.ll) 
Corollary A. 3 For every r/2 > fj, we have 

0(e;r/2,e,-e) < 1 for C G [fy, 7^2] (A. 12) 
for all0<e<l/2, and < e < 1. 

Proof of Corollary IA.31 A violation of (jA.12p implies the existence of 
rji € [?y,?/2) with 0(771) = 1 and 0'('?i) ^ Oi tlien (|A.11|I contradicts the condition 
0'(?72;'72,e,-e) = -e > -1. ■ 

Proof of Lemma IA.21 We fix 771 ^ 1 and define the variable y — C/771. 
Changing variables, equation ljA.8|) transforms into 

= ^ (0+^0) -10+^7/00 + 02-0, for 77>1 (A.13) 

vt \ y J 2 d 

0(l) = -i^7yi, 0(1) = 1, (A.14) 
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where ' = 4^. Define 

dy 

yo = sup{y > 1 : 4>{y) < and (f>{y) > 0} 

and note that yo > 1 if iji is large. On [1, j/o], < —Drji; therefore < 1 + l/(-D?7i), 
and consequently, on [l,?/o], 

1 , I fl 1\ , :, .„ . 8d 1 



1 



Similarly, 



Therefore, 



d+1 



< 



4 V 2 



101 



ml 



< 



1 / I 



4 V2 



if 



if 



m > 



d-2 dD 



d-2 



(d+l). 



1 •• 1 /I 1 



for large rji. Estimating \(j)\ by -D?7i, we find 



2 ' 



< 







1 


4 


V2 





[l,?/o], 



on [l,z/o], 



and since the right-hand side of this expression is negative for large 771 it follows that 
< on [l,?/o]; therefore yo may be redefined as 

yo = sup{y > 1 : (j){y) > 0}. 

It follows that on [l,yo]7 



< 



Dm 



1 1 

2 ~ 

1 1 

2 ^ 



d-2 
2d 

d-2 
2d 



ADm 



When < (f){y) < 1/2, this expression is bounded from above by —rif/64 for large 
rji. In terms of the original variable ^ we obtain 0'($) < —771/64, thus proving the 
lemma. H 



A. 3. Definition of p inTZ 

The general idea is to use rjo = fj as a reference point. In this way, owing to 
CoroUarv I A. 31 we can obtain the required estimates for p. It can happen, however, 
that the function (/)(^, ; 77, v, w) is not defined a,t ^ = fj. In such a situation, to define 
p, we introduce functions representing the intersection of 77, w, w) with the lines 
(j) — for T] < fj and (j) — 1 for rj > fj. Thus it is useful to define the following subsets 
of 7^ : 

T^i — {iVyVyw) S TZ: 4>{£,;ri,v,w) satisfies (i) in Lemma fA.ljk 

T^2 = {('^ji'jw) e TZ: r] , V , w) satisfies (ii) in Lemma FA. II} . with 

7^2a = 7^2 n {r; < ^} and 7^26 = 7^2 n {77 > 77}; 

TZs = {{ri,v,w) e TZ: (j>{£,;ri,v,w) satisfies (iii) in Lemma [A. Ij} -. 
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We treat the cases in turn. 

Case TZ^. Fix a point (j],v,w) G TZ^. We choose 770 = fj, vq = (t>(fj;ri,v,w) and 
= 4''{fj;r],v,w). Note that this choice is well defined: since — > C > as 
^ — > 00, there exists rje > fj such that (j){r]^;r],v,w) = e < 1/2, and —4>'{ri^]rj,v,w) — 
e < 1, with e ~ 1^. Then CoroUarv I A. 31 imphes that the solution (/){■; ri,v,w) can be 

continued to fj, even if ^ < 77. Setting p{rio, vq, wq) = Tyo^^e^''^/^, we find (c/. (|A.10|l ) 
Piv, V, w) = Ty'^+^e"'''/'' ejipl- £,(/){£,; f], (l){f]; rj, v, w), ((>' [f]; 77, v, w)) } . (A.15) 



fl 

The choice of 770 = ^7 also allows us to estimate the value of (p for S, > fj, which in turn 
permits us to control p for large 77, since the bound </>(^) < 1 for S, > fj implies an 
exponential decay for p as 77 — > cxd. 

Case TZi. Points in TZi are of the form (77, 1,0) and (7;, 0,0). We again choose 
770 = 77; substituting (j) = 1 and = into formula ljA.15|l gives 

p(77,l,0) =77'^+ie"^^^Tra^e"^, and p(?7, 0, 0) = T/'^+^e^'f . (A.16) 

Case TZ2a- Fix a point (77, v, w) G 'R-2a- Let rj* be given by Lemma lA.ll and define 
the function Lq: Ti2a such that LQ{ri, v, w) — min{7;*, 77}. Note that the function 

Lq is continuous and equals either the point 77* where 4'{r]* ; r] , v , w) vanishes or 77 if 
(j){fj;ri,v,w) > 0. To find p, we choose {rio,VQ,wo) = {rj* ,0, (/)' {r/* ,r],v,w)), and set 
p{r]o,vo,wo) = 77o+^e"''o/'*. This gives 

p(77, V, w) = 77^+1 exp{-r,V4 + /q}, (A.17) 

where 

I) 

/o = J {£.-, Loiv, V, w), (j)iLo{r], v, w);t], v, w), 0'(Lo(?7, w); 77, v, w)) d^. 

Lo{ri,v,w) 

Case TZ2b- Here it is convenient to define for any (r],v,w) G i?26 the function 
Ll:7^2fc ^R+, by 

max{77, max{^ € (0, 77) | (/)(^; 77, v, w) > 1}} if t; < 1, 

Li{t^,v,w) = { (A.18) 
min{^ e (77, 00) I (/)(^; r/, v, w) < 1} if t; > 1. 

The function Li is well defined for t; < 1 since if (/)(^; 77, v,w) = for some ^ € (77, 77) 
then (/) < 1 in (77,77) by Corollary I A. 31 and (p has to attain a local maximum in {fj,'i]), 
which is a contradiction with equation l)A.8(l . For t; > 1, Li is well-defined by Lemma 

ro 

Note that (f){Li{r],v,w);r],v,'w) < 1. The function Li is continuous and equals 
either 77^* where (j){rii:;r], v,w) = 1 or fj if <p{fj; 77, v, w) £ (0, 1). 

Now fix a point {r],v,w) £ Tl2b, choose 770 = Li{ri,v,w) and set p{r/o,VQ,wo) = 
^d+ig-(d-2),,^/4dg-^V2<i. Using (flToll . we find 

p(77, V, w) = 77^^+1 exp{-77V4 + f]l/2d - fj^ /2d + h} (A.19) 

where 

V 

j 2^(t){S,,Li{r],v,w),(j){Li{r],v,w);T],v,w),(t)'{Li{T],v,w),r],v,w)) d£_, 

Li{r]^v ,w) 
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and 

p(r], V, w) = 7?'^+! exp{-(d - 2)7/V4rf - f}^ /2d + I[} (A.20) 



where 

I[= J 2^[(l}{^,Li{i^,v,w),(t){Li{i^,v,w);r],v,w),(j}'{Li{r],v,w),r],v,w)) - l]d^. 
A. 4- Properties of p and ^ 

In the previous section, we have found a solution p of (|A.4|I . Here we show that this 
solution, together with the function $ given by l)A.5p . satisfies the properties required 
for the proof of Theorem 15.11 We start by stating a result which provides a lower 
bound for p in TZ2b- 

Lemma A. 4 Let M and M be the constants in estimates \44\l arid j^5| ), and let Li 
be defined as in iA.lS\) . Then there exists a large constant fjo such that the function 
G: [f]o,oo) — !■ given by 

G'(77) = max{ii(ry,a, -&) I 1 < a < A/ andO < < Af } for rj > fjo, 

satisfies G{r]) < Crj for some constant C = C{M) > 0. 

Proof. We take fjo large and we fix 77 > fjo. Using the continuity of Li, we have that 
G{ri) — Li{ri,d, —b) for some d S [l,Af], and b G [0,M]. Now we define the variable 
y = ^/i] ^ 1; the result is proved if we show that sup{y > 1 : 4>{y) > 1} < C{M). 
As in the proof of Lemma equation (|A.8|I transforms into 

0(1) = -bT] and (t>{l) = d. (A.22) 

Note that for ip > 1 we have (l){y) < for all y > 1, since (p{y) ~ implies that y can 
only be a maximum, a contradiction with equation ljA.8|l . 

We prove the claim in two steps. In the first step we consider the case 
d > d/2 - S > 1, where S = {d-\- l)d/ri^. Define yi = sup{2/ > 1 : ^(y) > d/2 -5}. We 
write (|A.21|I as 

^ -yA2{y)^ - Ai{y) for y > 1, (A.23) 

where 

'1 , 1 



Q^-[4>+— 4>]-7;4>+-vH + -4^ for y>l, (A.21) 



^i(2y)-r-0 and A2{y) 



d^ 2 ' y'^rf 



Since </>(•) G [d/2 — S, a] on [1, yi], A2 is non-negative and bounded by A2 := a/d. The 
function Ai is positive and bounded from below: 



d fd 



A,{y)>M:^{^-S -l--S]>0. 



Integrating equation ljA.23|) . we have 



y 



-r)^ I tA2(t) dt f -r)^ f tAiit) dt 

(t){y) = ~br]e 1 Ai{s)e = ds for y > 1. (A.24) 
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We observe that 

o } -V^JtA2it)dt 

Ai{s)e = ds>Ai^f{y;r)) for 1 < y < yi (A.25) 



1 

where f{y;ri) ~ if J e^^ '^^'■^ ds is a positive bounded function satisfying 

1 

yfiViV) 1/^2 as y ^ oo (the latter claim follows from considering the integrand 
close to s = y), and more precisely, 

y/(y; rj) > — =- for y > 2 and for sufhciently large rj. 

Therefore the primitive function 

F{y;v) = f{s;v)ds 

satisfies 

^^(2/;'7)> ^(logy-log2). (A.26) 

Integrating (|A.24p on [l,yi] and using ljA.25p . we obtain 

Hvi) < a-bri{yi - I) - AiF{yi]ri). 
To obtain a bound on yi, we use 0(yi) = d/2 — S and conclude 

AiF{yi;Tj) <d<M, 

from which it follows that yi < C(M) by (|A.26|) . 

For the second step, we replace rj by yirj in the rescaling above, by which we 
can assume that we are in the same situation: 0(1) = a, 0(1) = brj, but this time 
l<d<d/2-6. 

Similarly define y2 = sup{y > 1 : (j){y) > 1}. Since 1 < 0(-) < d/2 — S on [l,y2], 
the function A2{-) in ljA.23|) is negative, so that satisfies the differential inequality 

^0<-0' + 0< -2(0-1). (A.27) 

Let the function ijj solve 

= -2(V' - 1), with V(l) = a and V'(l) = 0. 

V 

The solution of this equation is V'(y) = 1 + (a — 1) cos{r]y/2{y — 1)), and note that 
■0(2/2) = 1 for y2 := n/{2r]y/2). From HA.27|) . 0(1+) < tp{l+)] if 0(y) = ipiy) for some 
y € (Ij y2), then by the comparison principle (which the operator ui-^ il / rj'^ -\- 2{u — 1) 
satisfies on intervals of length less than y2) we find > -0 on the interval [l,y], a 
contradiction with the previous remark. 

In conclusion we find that y2 < y2, thus proving the lemma. ■ 



We now derive estimates for p and in 7?. and TZ. 
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Lemma A. 5 The function p is continuous in TZ \ {i] = fj, v > 1}; for (j],v,w) G TZ, 
one finds 

P{t^, V, W) < T^d+l^-{d-2W/id_ (A 28) 

In addition, if {rj, v, w) G TZ, then 

P{V, V, w) > -^rf+^e-~^°''^ (A.29) 
for some constant Co — Cq{M) > 0. 

Proof. We start by proving (f08)l - lf09)l . Let = U n Ui for i = 1,2,3. If 
(77, V, w) € TZi then the estimates (IA.28|l - (jA.29|) follow by definition. If (77, v, w) G 'R-2a, 
then as (/) > on (77, LQ{r], v, w)) the integral in (|A.17|) is negative. This gives 

/3(?7,'f, w) < jy'^+^e"'''/'^ for (77, w, w) e 7^2a• 

Now for {ri,v,w) £ 7^2a, we have that [rj, Lo{t],v,w)] C [0,7^], v £ [0,M] and 
w £ [—M, 0]. Then the continuity of i/i on [77, io(^, w)] implies that 

i-iLoir], V, w), (l){Lo{r], v, w);t], v, w), (l)'{Lo{r], v, w); 77, v, w)) \co{[rj,Lo{v,v,w)]) < Co 

where Co = Co{M, M,fj). Using this bound to estimate /q in ljA.17|l . we find 

C(M)?7'^+^e"'''/'^ < p{ri, v, w) for (r/, v, w) £ 

with C(M) < 1, since we have integrated backwards. 

For any {ri,v,w) £ TZ2bi we use (|A.20(I and find the upper bound 



2d 



p(T],v,w)<Tf^\-'^^-^^'^"^^''e-'^"l 

This estimate follows from the negative sign of the integral in (|A.20|I . In fact, for 
^ £ (77, ^1(77, w, w)), we have 

f (l>{C,Liir],v,w),(j){Li{r],v,w);rj,v,w),(t>'{Li{rj,v,w);ri,v,w)) - 1 > if u > 1, 
\ (l>{CiLiir],v,w),(j){Li{r],v,w);rj,v,w),(t>'{Li{rj,v,w);r],v,w)) - 1 < if u < 1. 

Next for {r/,v,w) £ TZ2b, we find 

J p(ry, V, w) > 77''+ ^e"'''/'^ for v < 1, 

I p{v, V, w) > r^d+i^~{d-2)e/M^-ny2d^-ciM)n- „ > ^ 

where C{M) > 0. The estimate when v < 1 follows directly from l|ITT9)l . To obtain 
the estimate for p when w > 1, we use (|A.20|l . In fact, noting that (j) is non increasing 
in [77, Li{r], v, w)], we have that 

\(f) (•; Li(?7, V, w), 0(£i(7y, v, tw); 77, v, w), 0'(Li(7/, v, tw); 77, v, w)) |co([»?,Li(»).i;,tu)]) < M. 

Using this bound together with the estimate Li{ri,v,w) < C{M)r] (see Lemma lA.4|) . 
we find that /( in (|A.20(I satisfies — /( < C{M)rf' , which gives the derired estimate. 

To prove (|A.28|I for {ri,v,w) £ TZs, we examine two cases, if 7; < 77 then the 
estimate for Tl2a holds and for rj > fj the estimate for TZ2b holds. Finally, to obtain 
HA.29|) for (7;, V, w) £ TZs, we also check two cases, if 7; < ^ then the estimate for 7?,2a 
holds and for 7/ > 77 the estimate for TZ2b with 7; < 1 holds. 
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Claim, p is continuous in TZ \ {rj — fj , v > 1}. 

Before we prove this, note that TZ2 is an open set and TZi and TZ^ are cfosed. 

We first see that p is continuous within TZ2a and 7^2f): by continuity of Lq and Li. 
For the elements in TZi, the definition of p is as for 7?.2, therefore there is continuity 
of p between TZ2 and TZi. 

The dehcate part is to proof continuity between TZ^ and 7^2- Taking a 
sequence {r]n,Vn,Wn) G 7?,2, we associate a solution 4'ni-,'iln,Vn,Wn)- Suppose that 
iVn, Vn, Wn) w) £ 7?.3. Now if 7?, w, w) is the solution of (|A.8|I then ^ (/> 

in compact subsets of R+ . Therefore by Corollary IA.3L for n > no S N, we find 
4>n(jj) € (0, 1). Then {rin,Vn,Wn) S 7?.2 for n > no, have the same definition of p as 
for (77,1;, w) S 7^.3. Finally if w < 1 and rj — fj, then p is continuous. If 77 close enough 
to fj then we have that rjo = fj. So the computation of p uses the same formula, 
independent of the subset oiTZ to which [rj^v^w) belongs. H 

For <i> we deduce the following lemma, which implies (|56|l . 

Lemma A. 6 The Junction <i> is continuous inTZ\{ri — fj, v > 1} and if {rj, v, w) e R, 
then 

$(ry, 1;, ^) < |w;2 + ^d+l^-id-2W/4d 

and 

Hv,v,w)>-!^^-^-^^r,'+'e-('-'^^'/''. 

Proof. Follows directly from the definition (jA.5p of $ and uses the upper bound 
l|X28|) of p. U 



A. 5. Regularizing argument 

In the beginning of this appendix, we formally constructed a Lyapunov functional 
E{t) with $ and p satisfying (|A.3(I . In the previous section, we obtained a solution p 
of (|A.4p and $ given by (|A.2|I . Moreover these functions satisfy the properties found 
in Lemmas IA.5I and I A. 61 From these results we do not obtain enough regularity to 
derive (|A.3|1 . To do this, we introduce a introduce a regularization of $ using standard 
moUifiers and translation function to avoid the singularity of / at 77 = 0. See the details 
of the proof in [13 p. 102]. 



B. Appendix: Linear stability of blow-up profiles 

In this appendix, we study the linear stability of the blow-up profiles ipi and ip*, see 

Let i? be a solution of 127|I - H29|I and let ip he a solution of H19|l. The idea is to 
study the linearized equation for the difference $(77, t) := i?(?7,r) — ^{ri), i.e. 

d+l /I 1\ /I \ , , 

*r = %r, + + (j^ - 2 ) + b^'^" +2^-1] (B.l) 

Here, we have implicitly assumed that sufficiently close to blow-up only the linear 
terms play a role in describe the singularity formation. 
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For the stabihty analysis, let A > and consider a solution of HB.1|I of the form 
ipx{r))e^^. By (|RT1) . {iPx{v),>^) satisfies 

(■0a)w + ^^^(V'a)t, + Qiys - 'yCV-A),, + Qwr, + 2^ - 1 - Va- (B.2) 

For the analysis of boundary conditions we consider first ip — ipi. We note that at 
77 = we have either tp\ ^ 1 or tp\ ^ 1/rj'^. To have ipx bounded near 0, we impose 

[iPx)n{v)^^, as 77^0. (B.3) 
For large 77, we can either have tp\ ~ 7^^(2^+3)gi7 /4 qj. ^ rj'^^^^. We see that both 
behaviours diverge with 77, however the second asymptotic is bounded in terms of r 
and t as t T. Therefore to have a polinomial behaviour at infinity, we prescibe 

■0A(?7)e"'' 0, as 77 — > cx). (B.4) 

Now solving equation (|B.2(I together with (|B.3(I and ljB.4|l . we find a sequence solutions 
of HB.1|I given by {e^"'^Vn(»7)}neNu{o}i with Aq > Ai > . . . , where ipn ■= '0a„- If the 
blow-up time T > is chosen correctly in the definition of 77 and t, we can eliminate, 
see 201 1 the first mode (77 — 0) corresponding to change of blow-up and write 

B(77, r) = (^(77) + MV^'" + Oie^n- 

Therefore from the sign of Ai we obtain the linear stability of ip. 

In jnii Brenner et al. proved, using l|B.l(l . the following stability result for various 
blow-up profiles. 

Theorem B.l Every solution ip of ^4-7)) satisfying rjipfj/ip ^ 2 as 77 00 has an 
unstable mode corresponding to changing the blow-up time. Also, a blow-up profile 
with k intersections with the singular solution ips has at least k—1 additional unstable 
modes. 

In addition, the authors in found numerically that Ai < when ip — ipi and 
d > 2. In particular, they computed Ai = —0.272 ... for d = 3. This implies that (pi is 
linearly stable for d > 2. 

For (p = ip* , we can proceed as above and solve the eigenvalue problem for IjB.ip . 
Considering ljB.2|) with ip — ip* , we find that (■0a, A) satisfies 

{^>.)vv + - {^x)v + (1 - ^)V^A = 0, (B.5) 



Tj 2d 

with 1B.3|I and HB.4() . These boundary conditions are chosen by the same arguments 
for ip — ipi] however in the current case we have either ~ 7^d^(-*'^i)^'^~2g < 'jj^ 
or il)x ~ r;^^^^~^^ as 77 — > 00. Note that by changing rj by (—77) the equation remains 
invariant, so only solutions consisting on even powers are allowed. Then we construct 
a sequence of solutions of the form 

n 

il^niv) ^^A^jf' for any 77 = 0,1,2,3..., 

where the coefficients are given by Ai{2i{2i-\) + {d+l)2i) ^ A^^i{l-X-2i{d-2)/2d) 
for 7 = 1,2,... and Aq an arbitrary constant. This means that when (1 — A — 2(77. -|- 
l)(d — 2)/2d) = 0, we find an explicit polinomial solution of degree 277, where A is 
given by 

d-nid~2) 

A„ = . (B.6) 
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Consequently, we have obtained an explicit sequence of solution {(^n, A„)}„gfjy|o} 
for the eigenvalue problem (|B.5p . The eigenvalue Ag = 1 corresponds to the unstable 
mode of change of blow-up time and since Ai > for all d > 2, by ljB.6p . this means 
that if* is linearly unstable. 

References 

[I] J. Bebernes, and D. Eberly, A description of self-similar blow-up for dimension N > 3, Ann. 

Inst. H. Poincare, Anal. Nonlin. 5 (1988), 1-22. 
[2] M. D. Betterton, and M. P. Brenner, Collapsing bacterial cylinders, Phys. Rev. E 64 061904 
(2001). 

[3] P. Biler, Existence and nonexistence of solutions for a model of gravitational interaction of 

particles. III, Colloq. Math. 68 (1995), 229-239. 
[4] P. Biler, D. Hilhorst, and T. Nadzieja, Existence and nonexistence of solutions for a model of 

gravitational interaction of particles, II, Colloq. Math. 67 (1994), 297-308. 
[5] P. Biler, A. Krzywicki, and T. Nadzieja, Self-interaction of Brownian particles coupled with 

thermodynamic processes. Rep. Math. Phys. 42 (1998), 359-372. 
[6] P. Biler and T. Nadzieja, Existence and nonexistence of solutions for a model of gravitational 

interaction of particles, I, Colloq. Math. 66 (1994), 319-334. 
[7] P. Biler and T. Nadzieja, Growth and accretion of mass in an astrophysical model II, Applicationes 

Math. 23 (1995), 351-361. 
[8] M. P. Brenner, L. S. Levitov, and E. O. Budrene, Physical mechanisms for chemotactic pattern 

formation by bacteria Biophysical J. 74 (1998), 1677-1693. 
[9] M. P. Brenner, P. Constantin, L. P. Kadanoff, A. Schenkel, and S. C. Venkataramani, Blowup in 

the Chemotaxis equation. Technical Report preprint (1998). 
[10] M. P. Brenner, P. Constantin, L. P. Kadanoff, A. Schenkel, and S. C. Venkataramani, Diffusion, 

attraction, and collapse, Nonlinearity 12 (1999), 1071-1098. 

[II] C. J. van Duijn, I. A. Guerra, and M. A. Peletier, Global existence conditions for a non-local 

problem arising in statistical mechanics, submitted to Adv. Diff. Equations. 
[12] M. Escobedo, M. A. Herrero, and J. J. L. Velazquez, A nonlinear Fokker-Planck equation 

modellig the approach to thermal equilibrium in a homogeneous plasma, Trans. Amer. Math. 

Soc. 350(10) (1998), 3837-3901. 
[13] A. Friedman and B. McLeod, Blow-up of Positive Solutions of Semilinear Heat Equations, 

Indiana Univ. Math. J. 34(2) (1985), 425-447. 
[14] V. A. Galaktionov, On asymptotic self-similar behaviour for a quasilinear heat equation: single 

point blow up, SIAM J. Math. Anal. 26 (1995) 675-693. 
[15] Y. Giga and R. Kohn, Asymptotically self-similar blowup of semilinear heat equations. Comm. 

Pure Appl. Math. 38 (1985), 883-901. 
[16] Y. Giga and R. Kohn, Characterizing blowup using similarity variables, Indiana Univ. Math. J. 

36 (1987), 1-40. 

[17] I. A. Guerra, 2003 Stabilization and blow-up for some multidimensional nonlinear PDE's PhD 
Thesis TU/e 

[18] A. Haraux, 1991 Systemes Dynamiques Dissipatifs et Applications, Recherches en Mathemati- 

ques AppUquees, RMA 17, MASSON. 
[19] M. A. Herrero, E. Medina, and J. L. L. Velazquez, Finite-time aggregation into a single point in 

a reaction-diffusion system, Nonlinearity, 10 (1997), 1739-1754. 
[20] M. A. Herrero, E. Medina, and J. L. L. Velazquez, Self-similar blow-up for a reaction-diffusion 

system, J. Comput. Appl. Math. 97 (1998), 99-119. 
[21] W. Jager and S. Luckhaus, On explosions of solutions to a system of partial differential equations 

modeUing chemotaxis. Trans. Am. Math. Soc. 329 (1992), 819-24. 
[22] E. F. Keller and L. A. Segel, J. Theor. Biol 26 (1970), 399-415, 

[23] H. Matano, Convergence of solutions of one-dimensional parabolic equations, J. Math. Kyoto 

Univ. 18 (1978), 221-227. 
[24] J. Matos, Convergence of blow-up solutions of nonlinear heat equations in the supercritical case. 

Proc. Roy. Soc. Edinburgh Sect. A 129(6) (1999), 1197-1227. 
[25] J. Matos, Self-similar blow up patterns in supercritical semilinear heat equations, Commun. 

Appl. Anal. 5(4) (2001), 455-283. 
[26] A. A. Samarskii, V. A. Galaktionov, S. P. Kurdyumov, and A. P. Mikhailov, Blow-up in Problems 

for Quasilinear Parabolic Equations, Walter de Gruyter, Berlin, 1995. 



Self-similar blow-up 



28 



[27] J. J. L. Velazquez, Classification of singularities for blowing up solutions in higher dimensions, 

Trans. Amer. Math. Soc. 338(1) (1993), 441-464. 
[28] J. J. L. Velazquez, Stability of some mechanisms of chemotactic aggregation. SIAM J. Appl. 

Math. 62 (2002), no. 5, 1581-1633 
[29] G. Wolansky, On steady distributions of self-attracting clusters under friction and fluctuations, 

Arch. Rational Mech. Anal., 119 (1992), 355-391. 
[30] T. I. Zelenyak, Stabilization of solutions of boundary value problems for a second order parabolic 

equation with one space variable. Differ. Equat. 4 (1968), 17-22. 



